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ABSTRACT 

We estimate the stellar masses of disk galaxies with two independent methods: a photometrically 
self-consistent color-mass-to-light ratio relation (CMLR) from population synthesis models, and the 
Baryonic Tully-Fisher relation (BTFR) calibrated by gas rich galaxies. These two methods give 
consistent results. The CMLR correctly converts distinct Tully-Fisher relations in different bands 
into the same BTFR. The BTFR is consistent with Mb oc Vj over nearly six decades in mass, with no 
hint of a change in slope over that range. The intrinsic scatter in the BTFR is negligible, implying that 
the IMF of disk galaxies is effectively universal. The gas rich BTFR suggests an absolute calibration 
of the stellar mass scale that yields nearly constant mass-to-light ratios in the near-infrared (NIR): 
0.57 Mq/Lq in K s and 0.45 Mq/Lq at 3.6/im. There is only modest intrinsic scatter (~ 0.12 dex) 
about these typical values. There is no discernible variation with color or other properties: the NIR 
luminosity is a good tracer of stellar mass. 

Subject headings: galaxies: evolution - galaxies: fundamental parameters — galaxies: photometry 
— galaxies: kinematics and dynamics — galaxies: stellar content 


1. INTRODUCTION 

The stellar mass of a galaxy is one of its most fun¬ 
damental characteristics. Based on our knowledge of 
stellar evolution, we expect to be able to use measured 
galaxy luminosi t y and color to estimate s t ellar mass (e.g., 
Bell &; de Jond 20011: iBell et al.1 120031 IPortinari et al'l 


20041 iZibetti et all 120091 Unto fc Portinaril 120131) . How¬ 
ever, population synth esis models remain uncertain 
(1 Conrov fc Gunnl 1201011. and are not always self- 
consistent ( McGaugh fe Schombertl I2014H . 

Despite considerable progress, a variety of difficult is¬ 
sues persist. The star formation history of unresolved 
populations is difficult to uniquely constrain. The con¬ 
tribution to the integrated luminosity of bright but short¬ 
lived stars in late stages of evolution remains uncertain 
(e.g., AGB stars : |Maraston|j2005l iMarigo e t al.1 120081 
Kriek et all 120101 [Melbourne et al l 120121 : IZibetti et all 


20131) . The distribution of stellar metallicity as well 


as t he mean value appears t o play an important 
role dSchombert fc Rakosl 12009 : iConrov fe Gunnl 120101 
ISchombert fc McG augh 2 014aU bl). Stellar mass estimates 
can depend considerably on the treatment of these ef¬ 
fects. 

The largest persistent uncertainty is in the mass spec¬ 
trum of stars (the IMF). Most of the light is produced 
by massive stars while most of the integrated mass re¬ 
sides in low mass stars. Consequently, small differences 
in the IMF can have large effects on the stellar masss-to- 
light ratio, Y*. Such variations would l ead to substan¬ 
tial s catter in the Tully-Fisher relation dTullv fc Fished 
iM3), which not observed. This implies that t he galaxy- 
averaged IMF does not vary wildly dMcGaugh et al.l 
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l200fl lVerheiIenlf2?i?ViT) . 

Methods to estimate stellar mass that are indepen¬ 
dent of the IMF would be helpful. One approach to this 
problem is to measure the vertical velocity dispe r sion i n 
face-on disk galaxies (e.g., iBershadv et al.l [20101120111) . 
Another is provided by the Baryonic Tully-Fisher Rela¬ 
tion (BTFR). This is one of th e stro n gest correlations 
in extragalactic astronomy dMcGaughl 12005 :1. providing 
an empirical link between mass and observed rotation 
velocity. 

Re cent work on gas d o minate d g alaxies dBegum et al.1 
20081 I Swaters et all 120091 iStark et all 120091 
Trachternach et al.1 120091 ) has made it possible to 
obtain an absolute calibration of the mass scale of 
the BTFR t hat is largely independent of stellar mass 
estimates (lMcGaughll2012l) . When the gas mass exceeds 
the stellar mass, the systematic error in the stellar mass 
induced by the IMF, etc., becomes a minor contributor 
to the error budget (lMcGaughll2011l ). The assumptions 
necessary in population models cease to be the dominant 
factor in the calibration of the BTFR. 

Here we examine a sample of galaxies with high qual¬ 
ity, extended rotation curves from 21cm interferometry 
and photometry in the optical and near-infrared ($2|). 
We first construct the luminous Tully-Fisher relation 
for each band, and the corresponding stellar mass 
Tully-Fisher relation and BTFR using self-consistent 
stellar population mass-to-light estimators (Q. We 
then estimate the stellar mass of each galaxy using the 
BTFR independently calibrated by gas rich galaxies 
(0, and use this as a check on the IMF (0. We 
compare these results and discuss their consistency (or 
lack thereof) with other methods (0, and summarize 
in m 
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2. DATA 

The data used here are based on the sample of galaxies 
with extensive photometry in o ptica l and n ear-i nfr ared 
NIR) bands assembled by IMcGaugh fe Schombertl 
2014). All galaxies employed here have been observed 
at B , V, and [3.6]. Many, though not all, also have /, 
J, and JG-band data . Many SINGS (|Kennicutt, et, alJ 
120031: iDale et a ll I2005D galax ies observed as part of the 
THINGS ([Walter et, al.ll2008li HI survey are utilized. The 
samp le further inclu des many new Spitzer observations 
(ISchombert &; McGaughll2014bD . These data are used to 
estimate the luminosity a nd ste ll ar mass o f each gal axy . 

We select galaxies from lMcGaugh fc Schombert : 201 1: 
for which credible measurements of the outer, flat rota¬ 
tion velocity Vf are available from 21 cm interferome¬ 
try. This measure provides the best available estimate of 
the characteristic rotation velocity of disk galaxies in the 
sense th at it minimizes the scatter in the Tully-Fisher 
relation dVerheiienl 120011 : IMcGaughl I2005D . The more 
commonly employed line-width is less precise, in that 
it results i n a larger irreducible scatter in the BTFR 
dVerheiienl 120011 : IMcGaughl 1201 2t iZaritskv et, al.l I2014H . 
This difference is critical for estimating stellar masses 
from the BTFR, as uncertainty in the rotation velocity 
propagates strongly into the implied mass. 

We selected gala xies for_ observation with Spitzer 
(jSchombert fc McGaughil2014b ) to extend the dynamic 
range over which the BTFR could be constructed. Galax¬ 
ies with existing high quality HI rotation cur ves were se¬ 
lected to supplement the THINGS sample dWalter et alJ 
I2008T) at both ends of the BTFR. The resulting sam¬ 
ple extends to both greater and lower luminosity than 
present in THINGS. The aim was to cover as much range 
in mass as possible with evenly sampled bins. This is 
done to avoid the common selection effects afflicting mag¬ 
nitude limited samples, which typically have many L* 
galaxies but few or none with L < 10 9 Lq. The result is 
a sample of galaxies that probes the BTFR over a large 
dynamic range in a uniform way, free of the usual bias in 
which the region around L* is probed many times while 
either end of the relation is under-represented. 

The emphasis here is to uniformly sample as much of 
the BTFR as possible with the highest quality data avail¬ 
able. In this vein, we require that galaxies be sufficiently 
inclined that sin (i) corrections are modest (i > 45°; see 
discussion in [Stark et ahl l2009f ). This excludes a num¬ 
ber of THINGS galaxies with otherwise good data (e.g., 
NGC 6946). We also exclude two galaxies for which we 
do not entirely trust the data: NGC 2841 and NGC 2903. 
NGC 2 841 is a w ell-known outlier from Tully-Fisher re¬ 
lation ( Bottema et al.l 12002 ). Whether NGC 2841 is a 
true outlier, or some systematic error is to blame, is be¬ 
yond the scope of this paper. However, it is not simply 
at one extreme of the scatter; it sits well away from all 
other data. We therefore suspect some systematic error 
is to blame. In th e case of NGC 2903 , the [3.6] pho¬ 
tometry is suspect. Ide Blok et alJ ( 2008 ) note problems 
constructing a mass model, and the resulting [3.6] magni¬ 
tude makes NGC 2903 an outlier from both the (V— [3 .6], 
[3.61) color-magnitude relation (ISchombert fe McGaughl 
l2014al) and the [3.6] Tully-Fisher relation, even though it 
does not deviate from the Tully-Fisher relation in other- 
bands, including K s . We therefore infer that the [3.6] 


magnitude of NGC 2903 is not reliable, so we do not 
consider it further. 

The final sample contains 26 galaxies spanning the 
range 50 < Vf < 300 kms -1 , 3 x 10 7 < M* < 
3 x 10 11 M 0 , and 10 s < M g < 5 x 10 10 M 0 . It spans four- 
decades in luminosity, while magnitude lim ited Tully- 
Fisher samples typically span two (e.g., iPizaeno et alJ 
120071) . This large dynamic range is more important for- 
constraining the slope of the BTFR than is the sheer 
numbers of galaxies. 

Sample selection is driven by data quality: it is con¬ 
structed to be as free of systematic errors as possible. Of 
course, one can never completely exclude the presence 
of systematic errors. The two that concern us most are 
residual systematic errors is in the adopted distances and 
in the inclination determinations. 

Direct distances measurements vi a Ce ph eids or the 
TRGB are adopted when available dTullv et all 12009 ). 
When no direct measurement is available, we adopt a 
Hubble flow estimate consistent with independent lo¬ 
cal g alaxy data ( H 0 = 75 kms -1 Mnc -1 : fSorce et all 
120121) . Galaxies lacking direct measurements tend to be 
the most distant, luminous galaxies or the less studied 
dwarfs — i.e., those at the extremes of the relation. Di¬ 
viding the sample into those objects with and without 
direct distance determinations reveals no statistically sig¬ 
nificant difference: we recover the local distance scale. If 
instead Hq = 67 kms -1 Mpc -1 , the primary systematic 
effect would be to increase our mass estimates by 0.1 dex 
with a corresponding translation of the BTFR zero point. 

The effects of syste matic inclination errors are dis¬ 
cussed at length bv IMcGaughl (|2012D . These most com¬ 
monly manifest as outliers to the low velocity side of 
the BTFR. Oval distortions can make galaxies look less 
face-on than they are, so the inclination correction may 
be underestimated. This is mitigated here by the use 
of velocity fields, which give an additional constraint on 
the inclination. As we will see, the scatter around the 
BTFR relation is small with no obvious asymmetry, so 
systematic inclination errors appear to be small in this 
sample. 

Fundamental d ata lik e dis ta nces, magnit udes, and col¬ 
ors are given in iMcGaugh fc Schombertl ( 2014 ). Data 
specific to our present analysis are given in Table [I] Col¬ 
umn 1 gives the name of the galaxy. Column 2 gives the 
rotation velocity Vf. Columns 3 and 4 give the atomic 
and molecular gas masses, respecively. Luminosities in 
each of the V, I , and [3.6 ] bands are given in columns 
5 7. As in IMcGaugh fc Schombertl ( 2014 ), we adopt 

solar absolute magnitudes of 4.83, 4.08, and 3.24 in the 
V, I, and [3.6] bands. References to the sources of the 
data are given in column 8. 

Rotation curves generally experience an extended ra¬ 
dial range over which the rotation velocity is approxi¬ 
mately constant. This flat velocity, Vf, is an obvious 
measure of the characteristic global rotation speed of a 
galaxy. We have measured or remeasured anew Vf for 
all galaxies considered here by taking the mean of data 
along the outer, flat portion of the rotation curve. While 
the outer parts of rotation curves are very nearly flat, 
they are not perfectly so. We implemented several algo¬ 
rithms for deciding what constitutes the flat portion of 
the rotation curve: a subjective estimate by eye, a defini- 








































































Weighing Galaxy Disks 


3 


Table 1 

Rotation Velocities and Gas Masses 


Galaxy 

v f 

(kms 1 ) 

M{HI) M(H 2 ) 
(10 9 M 0 ) 

L v 

Li 

(to 9 L 0 ) 

^[3.6] 

Refs. 

DDO 154 

53.3 ±4.7 

0.0925 

0.00618 

0.0516 

0.0294 

0.0725 

1,2,3,4,5 

D631-7 

52.9 ±5.0 

0.147 

0.0127 

0.0540 

0.0448 

0.0975 

1,2,3,6 

DDO 168 

53.4 ±2.5 

0.317 

0.00279 

0.166 


0.187 

1,2,3,7 

D500-2 

68.1 ±6.6 

0.871 

0.00790 

0.304 

0.226 

0.395 

1,2,6 

NGC 2366 

59.7 ± 10.3 

0.600 

0.123 

0.456 

0.369 

0.716 

1,5,8 

IC 2574 

76.7 ±4.5 

1.41 

0.0937 

1.04 

0.956 

2.33 

1,5,8 

F563-1 

111. ± 10. 

3.50 

0.0900 

1.13 

1.25 

2.84 

1,2,3,7 

NGC 2976 

86.0 ±3.5 

0.134 

0.121 

1.13 

1.05 

3.19 

1,4,5 

F568-V1 

124. ± 10. 

3.93 

0.134 

2.36 

2.25 

4.20 

1,2,3,7 

NGC 1003 

113.5 ± 1.9 

5.31 

0.243 

2.76 


5.91 

1,4,5 

F568-1 

116. ± 10. 

4.02 

0.141 

3.10 

2.98 

7.10 

1,2,3,9 

NGC 7793 

110.4 ±4.4 

0.761 

0.282 

2.99 

1.80 

7.55 

1,4,5 

UGC 128 

131. ± 10. 

6.51 

0.420 

4.49 

4.20 

11.3 

1,2,3,10,11 

NGC 2403 

135.8 ± 4.0 

2.52 

0.522 

3.86 

5.00 

12.1 

1,4,5,10.11 

NGC 925 

113.6 ±6.3 

3.99 

0.826 

8.35 

8.35 

15.8 

1,4,5 

NGC 3198 

148.9 ± 4.2 

0.134 

1.30 

12.4 

15.8 

31.2 

1,4,5 

NGC 3621 

152.3 ±3.2 

6.98 

7.06 

6.74 

7.11 

33.1 

1,4,5 

NGC 3521 

191.9 ± 17.6 

4.48 

1.65 

15.6 

23.2 

91.1 

1,4,5 

NGC 3031 

199.0 ± 13.8 

6.97 

1.58 

23.0 

38.7 

106. 

1,4,5 

NGC 5055 

181.1 ± 11.9 

7.21 

2.36 

26.3 

40.2 

135. 

1,4,5 

NGC 2998 

211.7± 2.7 

27.7 

7.27 

74.8 


160. 

1,7 

NGC 6674 

240.5 ±4.2 

38.2 

9.30 

63.2 


230. 

1,7 

NGC 7331 

245.8 ± 8.4 

9.38 

4.30 

38.3 

67.1 

259. 

1,4,5 

NGC 801 

218.8 ±3.5 

36.7 

12.0 

71.5 


268. 

1,7 

NGC 5533 

239.9 ±5.0 

23.5 

13.9 

62.7 


303. 

1,7 

UGC 2885 

298. ± 10. 

34.5 

15.7 

179. 


468. 

1,7 


Refe renc es. - 1. McGaugh V Schombert (2014). 2. Schombcrt & Mc Gaug h 

J2014bl). 3. ISch ombert et al.l d201U). 47 de Blok et a l. (2008). 5. Leroy ct al. (2008). 
6. ITrachternach et al.l (|2009l). 7. McGaug il (120051). 8. lOh et a 1 ([20081). 9 . Ide Blok et all 


JT996T) . 10. Ide Blok fc McGaugB (fl99(D . 


Note. — Galaxie s are listed 
I McGaugh fc Schombert! ([2014D . 


11. IVerheiien fc de Biol 
order of increasing [3.6] 


luminosity as in 


tion by slope such that (I log V /0 log R < 0 .1 (|McGaughl 
120051 : iSwaters et al.ll2009t IStark et ahll2009D , and an out¬ 
lier rejection routine that iteratively excludes points from 
the rising or falling portion of the rotation curve. These 
methods yield consistent results within the errors. In¬ 
deed, the relative variation around Vf along a rotation 
curve, while clearly perceptible, is typically smaller than 
the absolute uncertainty in Vf. AV = A R{dV/dR) < 
ay. Hence the mildly non-zero slopes dV/dR that are 
observed do not preclude accurate measurement of Vf. A 
greater concern is the occasional warp that causes a ra¬ 
dial variation in the inclination within a galaxy. Indeed, 
this might be one issue in NGC 2841. The galaxy for 
which this concern is most import ant in the present sa m- 
ple is DDO 154 (see discussion in Ide Blok et al.ll2008li . 

The gas mass of each galaxy is adopted from the cor¬ 
responding tabulated reference in Table |T] The HI mass 
is obtained from the 21 cm flux in the usual fashion: 
M(HI) = 2.36 x 10 5 D 2 Fhi- The molecular gas mass has 
been e stimated for the THINGS galaxies bv ILerov et all 
(j2008l) . who discuss in detail the conversion from ob¬ 
served CO flux to H 2 mass. 

Many galaxies lack CO detections. In these cases, 
we exploit t he near-const a nt effi ciency of star formation 
observed by iLerov et~ahl (2 008 ) to estimate the molec- 
ular gas mass from th e observed star formation rate 
(ISchombert et al1l201 If) . This method implicitly corrects 
for the metallicity dependence of the CO-to-H 2 conver¬ 
sion factor, albeit by assuming a constant star formation 
efficiency. For galaxies lacking CO detections, we esti¬ 


mate the molecular gas mass as 

log[Af(ff 2 )] = log(SFR) + 9.15. (1) 

Here, the molecular gas mass is in solar masses and the 
star formation rate SFR is in M 0 yr -1 . These estimates 
are rather u ncertain (compare to other estimators dis¬ 
cussed by iMcGaughl 120121) , but the molecular gas mass 
is usually a small portion of the mass budget, and makes 
little difference to the result. 

The luminosities given in Table Q] are computed di¬ 
rec tly from t he distances, magnitudes, and colors given 
in iMcGaugh fc Schombert! (1201411. The [3.6] lumi- 
nosities originate either with Ide Blok et~ahl (|2008fl or 
ISchombert fe McGaughl (|2014bT ). with other bands as 
cited in the table. The formal uncertainties in the mea¬ 
sured magnitudes are small, of order 0.02 mag. The dom¬ 
inant uncertainty in the luminosity and the gas mass is 
the distance. We treat the uncertainties in the luminosi¬ 
ties as negligibly small, and consider the random distance 
uncertainties as a contributor to the scatter in the Tully- 
Fisher relations we construct. 

3. DIRECT FITS 

3.1. Types of Tully-Fisher Relations 

The term “Tully-Fisher Relation” has been used 
to mean a variety of subtly different things. Most 
commonly, the Tully-Fisher Relation is posed as a 
relation between ab solute magnitude and line-width 
(jTullv fe Fisherlll977lf and employed as a distance indi¬ 
cator (e.g., iSorce et al.ll2012f) . Here we construct several 
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Figure 1. Tully-Fisher relations determined from V (left), I (middle), and [3.6] (right) data for the galaxies in Table [T] The top row 
shows the luminous Tully-Fisher relation (LTFR), the middle row the stellar mass Tully-Fisher relation (STFR), and the bottom row the 
baryonic Tully-Fisher relation (BTFR). In all cases the luminosity or mass is plotted against the flat, outer vel ocity . Fits to eac h relation 
are shown as dashed lines and reported in Table(2] For reference, the BTFR, previously fit to gas rich galaxies bv lMcGaughl (|201 2i ) is shown 
as a light solid line in all panels. 


distinct flavors of Tully-Fisher relations from the assem¬ 
bled data in hopes of identifying the underlying physical 
basis of the relation. 

We start with the relation between luminosity and 
asymptotic rotation velocity in each band. This is anal¬ 
ogous to the common Tully-Fisher relation, but differs 
from it in that Vf from extended, resolved rotation curves 
is used rather than line-widths. While similar, these are 
not identical measures of the characteristic rotation ve¬ 
locity of a galaxy. The primary consequence is a system¬ 
atic difference in the slopes of the fitted relation. Line- 
widths are sensitive to the presence of inner humps in ro¬ 
tation curves, which are prominent in lu minous galaxies 
but largely absent in faint ones (see, e.g. lVerheiienll200TI : 
iNoordermeer fe Verheiienll2007H . This causes the L-W 
relation to be shallower than the L-Vf relation, which 


we will refer to as the LTFR. 

Next, we construct the stellar mass Tully-Fisher rela¬ 
tion (STFR). We employ stellar population models to 
estimate the stellar mass of each galaxy based on their 
observed luminositi es and the co l or-mas-to -light ratio re¬ 
lation (CMLR) of iMcGaugh fc Schombertl ( 2014 b. We 
then construct the STFR implied independently by the 
luminosity observed in each band pass. If the our semi- 
empirical CMLR are correct, then each band’s LTFR 
should transform into the same STFR. 

Finally, we construct the baryonic Tully-Fisher Rela¬ 
tion (BTFR). This is a relation between rota t ion veloc¬ 
ity a n d baryonic mass (IMcGaugh et, al.l [20001: iVerheiienl 
20011: iBell fc de Jond 120011 : iPfenniger fe Revaz! 120051 
McGaughl I2005H . where now the mass includes all ob¬ 
served forms of baryonic mass, stars and gas: Mb = 
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Table 2 

Tully-Fisher Relations and their Scatter 


Band 

i 

X 

log A 

^log I 

Luminous Tully-Fisher Relation 


V 

L v 

4.27 ±0.19 

0.68 ±0.41 

0.18 

I 

Li 

4.77 ±0.30 

-0.37 ± 0.62 

0.19 

3.6 

^[3.6] 

4.93 ±0.17 

-0.28 ±0.36 

0.15 

Stellar Mass Tully-Fisher Relation 


V 

M * 

4.70 ±0.22 

-0.20 ±0.48 

0.21 

I 

M* 

5.01 ±0.31 

-0.84 ±0.64 

0.20 

3.6 

M* 

4.93 ±0.17 

-0.61 ± 0.36 

0.15 

Baryonic Tully-Fisher Relation 


V 

M b 

3.92 ±0.18 

1.81 ± 0.39 

0.17 

I 

M b 

4.02 ±0.24 

1.56 ±0.51 

0.15 

3.6 

M b 

4.09 ±0.15 

1.49 ±0.32 

0.13 

Gas-only Tully-Fisher Relation 



M g 

3.24 ±0.28 

2.89 ±0.59 

0.29 


Note. — Fits to the data in Table [T] of the form 
log i = log A ± x log Vf. 


M* + M g . Again, if all is well, the same relation should 
be obtained independently from each band. We can also 
compare the result to that obtained independently from 
fitting gas rich galaxies, which is very ne arly i ndep endent 
of the stellar mass estimate (lMcGaughll201lL 2012 ;. The 
LTFR, STFR, and BTFR are shown in Fig. [T] 

The use of Vf as the velocity measure rather than 
line-widths leads to systematically steeper slo pes for 
all flavors of the relation. For the BTFR, [McGaughl 
(120121) found a slope x = 3.41 ± 0.08 when line-widths 
are employed as the velocity measure, c onsistent with 
the 3. 5 ± 0.2 reported independently by Za ritskv et al.1 
(|2014l) . In contrast, the slope of the BTFR is very close 
to 4 w hen Vf is employe d as the measure of rotation ve¬ 
locity (IStark et al.ll2009D . It also consis tently yields rela- 
tions with lower scatter. For example, [McGaugh (120121 ) 
estimates an irreducible scatter of ~ 0.2 dex in mass in 
the BTFR when line-widths are employed, but negligible 
intrinsic scatter with Vf. 

3.2. Fitting 

For the various possible ordinates of the Tully-Fisher 
relation, we fit relations of the form 

log(^) = log A + xlog(V f ). (2) 

Here A is the normalization of the relation (log A being 
the intercept in the log-log plane) and x is the slope. The 
symbol i represent any of the luminosity, stellar mass, 
gas mass, or baryonic mass, while Vf is of course the 
rotation velocity measured in the outer portion of the 
resolved rotation curve. 

To fit the data here, we use the maximum likelihood 
method of I Weiner et al.l (l2006i ). This allows for a finite 
intrinsic scatter in the ordinate in addition to random 
errors in both dimensions. The uncertainty that we es¬ 
timate in the abscissa Vf is given in Table [I] The error 
in luminosity is dominated by uncertainty in distance 
rather than in apparent magnitudes. The uncertainty 
in the distance to each galaxy is heterogenous and hard 
to quantify with confidence. We therefore treat it as a 
contributor to the scatter. The best-fit scatter found in 


this way is a combination of random errors and intrinsic 
scatter, placing an upper limit on the latter. The results 
of the fits are reported in Table [2] 

The slope of the LTFR becomes progressively steeper 
as we go to progressively red der bands (Table [p, consis¬ 
tent with previous findings dTullv et al.l 1 19981 IVerheiienl 
l200H) . As one might expect, the [3.6] LTFR has less scat¬ 
ter and better determined fit coefficients than that in V 
or I. Nevertheless, we caution against over-interpreting 
the specifics of these LTFR fits given the high gas content 
of the slower rotators. 

The LTFR found here are steeper than many published 
relations. This is expected, as the particular value of the 
fitted LTFR slope depends on the dynamic range of the 
data. As one goes to lower velocities, one finds more 
gas rich galaxies. The luminosity progressively underes¬ 
timates the mass, so these galaxies fall too low, steepen¬ 
ing the fitted slope. The scatter in luminos ity at a given 
velocity also goes up dMatthews et al.lll998D . The precise 
value of the fitted slope therefore depends on both the 
dynamic range of the sample, and the happenstance of 
the gas fractions of the galaxies therein. Magnitude lim¬ 
ited samples will be biased to lower slopes, as they will 
preferentially include brighter, lower gas fraction galaxies 
at a given Vf. This can be a strong effect, since it oc¬ 
curs at one end of the relation where relatively few points 
can have considerable leverage on the best-fit slope. In 
our sample, exclusion of the three slowest rotators (with 
Vf « 50 kins -1 ) flattens the fitted slope perceptibly (if 
not significantly). 

An important consequence of sample dependence is 
that there probably is no such thing as a universally cor¬ 
rect LTFR. The “right” answer is unique to each sample. 
This effect is difficult to discern in most samples, which 
tend to lack galaxies with Vf < 90 kms -1 where the 
effect becomes pronounced. This may be enough to ex¬ 
plain subtle differences in published calibrations. It may 
also introduce systematic errors in distance determina¬ 
tions, as one is sampling different parts of the relation at 
different distances. Fortunately, this effect will be small 
as long as only bright galaxies are utilized in both cali¬ 
bration and distance determination. 

3.3. Stellar Mass 

Stellar mass is a more fundamental physical quantity 
than luminosity. The LTFR can and does vary from 
band to band. These all presumably originate from the 
same STFR, with differences in the LTFR reflecting sys¬ 
tematic variations in the stellar populations composing 
fast and slow rotators. From the color-magnitude rela¬ 
tion, one would expect progressively steeper LTFR as 
one goes from blue to red bands, but these should all 
originate from one STFR. 

The stellar mass is estimated from M* = T*L for each 
band with the mass-to-light ratio Y* estimated by popu¬ 
lation synthesis models. Many models currently in wide 
use do not pass a straightforward self consistency test, 
in that stellar masses estimated from the luminosity in 
different pass ban ds for the same g alaxy return wildly 
different results ([McGaugh fc Schombertll2014D. By com¬ 
paring data over a broad range. [McGaugh fe Schomberfi 
(|2014l) constructed empirically self-consistent color-T* 
relations. 

We employ the CMLR of IMcGaugh fe Schomberfi 
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(l2014f) to estimate T* with B — V colors, which we found 
to be the most sensitive to variations in stellar mass-to- 
light ratio among t he a vailable colors. To be specific, we 
adopt the model of lBell et al.l 1 2003 1 modified to give self 
consistent photo metric r esults as quantified in Table 7 of 
IMcGaueh &; Schombertl (120141) : 

log = —0.628 + 1.305(£ - V) 
log Y* = —0.275 + 0.612(5 — V) (3) 

tI 3 ' 61 = 0.47 

As expected, the dependence of Y* on color becomes 
weaker as one goes to redder bands. It effectively disap¬ 
pears in the NIR, so we neglect the meaninglessly tiny 
color slope (—0.007) at [3.6]. 

We make this pa rticular choice of model because the 
model of lBell et al.l (j2003f l required the least modification 
of the various models considered. By construction, the 
U-b and Y*-indicator is identical to the original one of 
iBell et all (]2003l ). and is subject to the usual normal¬ 
ization uncertainty due to the IMF. Both the /-band 
and [3.6] Y* are modified to obtain photometric self- 
consistency. Requiring self-consistent stellar masses im¬ 
plies 0.4 < Y* 3 ’ 6 ^ < 0.5 for all models, including those 
that initially predicted rather different values. A con¬ 
stant NIR Y* is therefore a good approximation. 

The dependence of the NIR T* on color in eq. [3] is 
considerably weaker than the already weak relation an¬ 
ticipated in many models. This is required for self- 
consistency: the photometric data are not consistent 
withJN4R mass-to -ligh t ratios that depend much on color 
(see lMcGaugh fe Schombertll2014D . While there can cer- 

tainly be some color dependence, the variation of y 1 3 ' 6 ^ 
with color appears to be rather less than the intrinsic 

scatter in Yi 3 6 ^ from galaxy to galaxy, which itself is 
also fairly small. 

In practice, the weak relation between NIR Y* and 
color means that NIR luminosity is a good proxy for 
stellar mass. One might hope to do better by fitting the 
entire SED of a galaxy, and this no doubt matters for ex¬ 
treme objects like those dominated by the light of very 
young stellar populations. However, for most galaxies, 
there is no value added from this procedure. Most of 
the information about Y* is contained in the first color 
(B — V here), which is already negligible in the NIR. 
Further corrections from other colors carry progressively 
less information. In the case of [3.6] and the AVband, at¬ 
tempts at SED fitting are more likely to add noise than 
improve the estimate of Y*. Worse, they are prone to 
induce systematic errors from imperfections in the pop¬ 
ulation models. For example, many models assume a 
single metallicity for the stellar population. This ap¬ 
pears to be one driver for the predicted dependence of 
the NIR mass-to- light rati o o n co lo r bein g stronger than 
observed dMcGaugh fc Schombertl |20 141) . Models with 
realistic metallicity distributions h ave a much weaker de¬ 
pendence of NIR T* on color (|Schombert fe McGaughl 
I2014al) . 

If the stellar masses-to-light ratios estimated with eq. [3] 
are correct, the Y* for each band should succeed in con¬ 
verting the different LTFR observed for each independent 
band into the same STFR. This procedure is indeed suc¬ 


cessful. The STFR in Table [2] are equivalent within the 
errors. 

Since the adopted NIR mass-to-light ratio is constant, 
the [3.6] STFR is identical to the [3.6] LTFR except for 
the change in normalization due to the non-unity value 
of the mass-to-light ratio. Consequently, it is the slope 
of the V and / LTFR that change to result in a STFR 
that matches that from [3.6]. This goes in the expected 
sense from the color-magnitude relation: fainter galaxies 
are bluer and have lower Y* on average. Use of incorrect 
CMLR would lead to discordant STFR in the various 
bands. That we find consistent SFTR provides indepen¬ 
dent confirmation t hat t he CML R in eq. [3] determined 
photometrically by iMcGaugh fc Schombertl (I2014D are 
very nearly correct. 

The conversion from the LTFR to the STFR can only 
work to the extent that B — V is predictive of Y*. That 
it does work suggests that B — V is an adequate T* 
indicator. It certainly is not perfect: one expects a fair 
amount of scatter about the mean CMLR. The expected 
consequence is an increase in the scatter in the STFR 
relative to the LTFR. Even though the mean slope of the 
CMLR succeeds in converting the LTFR into the STFR, 
the scatter in the latter increases due to misestimates of 
the mass-to-light ratios of individual galaxies. That is, 
the scatter in the STFR should include whatever scatter 
is already in the LTFR, plus an extra component for 
the scatter in the CMLR. This effect does appear to be 
present in the data, as seen by the amount of scatter in 
Table [2] 

From a population perspective, one expects a greater 
amount of scatter in Y* at bluer wavelengths. The scat¬ 
ter fit in Table [2] is a combination of intrinsic scatter and 
that from random errors like those in distance. While 
there is considerable uncertainty in these uncertainties, 
we estimate that random errors contribute ~ 0.08 dex 
to the scatter in the ordinate of the various Tully-Fisher 
relations. Subtracting this in quadrature from the fitted 
value for each band gives an estimate of the intrinsic scat¬ 
ter in Y*: 0. 16 dex in V, 0.13 in /, and 0.12 in [0.36]. For 
comparison, IBell fc de Jond 1 2001 ]) anticipate 0.1 dex of 
scatter in K from variations in th e star formation history 
alone, and iPortinari et ahl (I2004D anticipate 0.15 dex. 

It thus appears that the observed scatter is entirely 
explained by random errors plus the expected scatter in 
mass-to-light ratios. The scatter in Y* is consistent with 
that expected from the inevitable variations in star for¬ 
mation histories. This leaves essentially no room for vari¬ 
ations in the IMF from galaxy to galaxy. Apparently, the 
IMF in disk galaxies is universal, at least averaged over 
the many star formation events that build a galactic disk 
over a Hubble time. 

There is also precious little room for intrinsic scatter 
in whatever physics underlies the Tully-Fisher relation. 
Given the small size of the sample, we should not over¬ 
interpret the precise values of the scatter. However, it 
is not obvious that larger samples will inevitably lead 
to larger scatter, as they bring with them the larger 
risk of un acco unted systematic errors (see discussion in 
lMcCaughll2012i ). Bear in mind that in order to increase 
the intrinsic scatter of this sample, we must have overes¬ 
timated random errors. Even in the limit of zero exper¬ 
imental error and no scatter in Y*, the intrinsic scatter 
is still small (< 0.15 dex; Table [2]). 
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Figure 2. The gas-only Tully-Fisher relation. Gas mass is plotted against rotation velocity without any reference to stellar luminosity. 
The data are plotted and fit (dashed line) in the left panel (see Table [2] for coefficients). The fit looks reasonable, though the scatter is 
larger than in any of the LTFR. The right panel shows the same data, but with points color coded by gas fraction. To estimate the ratio 
M*/M g , we use M* = 0A7L^ gi. Red points have M*/M g > 2, orange points 1 < M*/M g < 2, green points ^ < M*/M g < 1, and blue 
points M*/M g < i. The data clearly segregate by gas fraction. Gas rich galaxies follow the BTFR (solid grey line) while gas poor galaxies 
fall well below it. For comparison, the dotted lines parallels the BTFR but is shifted by a factor of five. This is the amplitude of the shift 
one would expect for a typical Milky Way-like spiral with a ~ 20% gas fraction. The segregation by gas fraction renders meaningless the 
fit in the left panel. 


3.4. Gas Mass 

The atomic and molecular gas masses given in Table [1] 
are used to compute the total gas mass 

M g = 1.33 [M(HI) + M(H 2 )\. (4) 

The factor 1.33 accounts for the mass in helium. One 
can debate the appropriate value of the second digit past 
the decimal, but not at a level that mat t ers h ere. We 
adopt 1.33 for consistency with iMcGaughl ( 2012 ). whose 
BTFR calibration will be utilized in ]J4] The contribu¬ 
tion of other gas phases is generally negligible, at least 
within the visible disk where all relevant quantities are 
measured. 

Luminosity is the original ordinate of the Tully-Fisher 
relation. Stars are the dominant baryonic mass compo¬ 
nent of bright galaxies. We therefore do not expect there 
to be a meaningful gas-only Tully-Fisher relation. For 
completeness, we do construct one. It is shown in Fig. [2] 

At first glance, the relation in Fig. [2] appears perfectly 
reasonable. It has considerably greater scatter than any 
of the STFR, but it is certainly possible to fit a line 
through the data (left panel of Fig. [2j. This is not par¬ 
ticularly meaningful, as the scatter here is real. Galaxies 
of the same Vf can differ in gas mass by an order of 
magnitude. 

Galaxies segregate by gas fraction (right panel of 

Fig. [2]). To estimate the gas fraction, we use xl 3 '^ = 0.47 
from eq. [3] to estimate the stellar mass, and bin the data 
by the ratio M*/M s . Galaxies with M* < M g hover near 
the BTFR previously fit to gas rich galaxies (IMcGaughl 
IMD- Since stellar mass is ignored here, these should 
all fall slightly below the BTFR line. Most do, though 
a few fall a tad too high, albeit within the errors. In 
contrast, galaxies with M* > M g all fall well below the 
line. The obvious inference is that we are missing an 
important component of the mass in these galaxies: the 
stars. Indeed, the amount of the shift is what we expect 
for the typical gas fractions of bright spirals. 


Just as we should not ignore the gas in faint galaxies, 
we should not ignore the stars in bright ones. Neither 
the LTFR nor the gas-only Tully-Fisher relation are fun¬ 
damental. Instead, they must stem from a more basic 
relation that does not distinguish between mass in the 
form of stars or gas. 

3.5. Baryonic Mass 

In order to construct the BTFR, we take the sum of 
stellar and gas mass to obtain the baryonic mass: 

M b = M* + M g . (5) 

The stellar mass is computed with eq. [3] and the gas mass 
is from eq. [H The BTFR resulting from the stellar mass 
as computed from each of V, /, and [3.6] is shown in the 
bottom row of Fig. [TJ 

If all has gone well, this BTFR should be consistent 
from band to band. This should necessarily follow if the 
STFR has been computed correctly for each band. In¬ 
deed, this is the case: the relations computed separately 
for each band are mutually consistent. Both the slopes 
and the intercepts are indistinguishable within the er¬ 
rors. This convergence is readily seen in Fig. [T] These 
relations are also consistent with the BTFR previously 
fit to indep endent d at a, for gas rich, predominantly slow 
rotators bv IMcGaughl (20121 ). 

Another encouraging point is that scatter in the BTFR 
is less than that in either the LTFR or the STFR. Includ¬ 
ing the gas mass as well as the stellar mass has made a 
better relation in this sense. Indeed, the relation is near 
perfect in the sense that their appears to be very little in¬ 
trinsic scatter. In 33.31 we estimated the intrinsic scatter 
in stellar mass-to-light ratios for this sample to be 0.16 
dex in V, 0.13 in J, and 0.12 in [0.36]. The scatter in 
the BTFR is 0.17 dex in V, 0.15 in J, and 0.13 in [0.36] 
(Tabled. These are effectively identical given the uncer¬ 
tainties, leaving precious little room for intrinsic scatter 
in the relation itself. 
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The fundamental, underlying physical relation appears 
to be one between baryonic mass and asymptotic rotation 
velocity. The data are consistent with a single power law 
(a single linear slope in log-log space) over four decades in 
baryonic mass. There is no evidence for a bend or break 
in the relatio n as so metimes anticipated by models (e.g., 
iTruiillo-Gomez et al.ll201lf) . 

The luminosity of the traditional LTFR is merely a 
proxy for baryonic mass. This works well for bright spi¬ 
rals that are gas poor, but starts to break down at lower 
velocity where galaxies tend to be more gas rich. As a 
consequence, the LTFR suffers increased scatter at low 
rotation velocity, and may appear to bend. This sim¬ 
ply illustrates the shortcomings of luminosity as a mass 
indicator, not any break in the underlying relation. 

Since luminosity is an imperfect indicator of baryonic 
mass, calibrations of the LTFR are inevitably inaccurate 
at some level. Bright spirals do contain gas, and the gas 
fraction varies from galaxy to galaxy. The LTFR will 
inevitably deviate from the underlying BTFR in a way 
that depends both on the band pass and the gas fractions 
of the sample of calibration galaxies. Fortunately this 
appears to be a minor effect for bright galaxies. 

The small intrinsic scatter in the BTFR suggests that 
there is no large reservoir of missing baryons remain¬ 
ing to be discovered in the disks of spiral galaxies 
(lMcGaueh et al.1 feOOOH . It is obvious if we miss either 
the stars or the gas (Figures [T] and [2]). If we were still 
missing a comparable number of baryons, then the sum 
in eq. [5] would be incomplete. This would show up as a 
systematic deviation from the BTFR, which is not appar¬ 
ent. One can always have some mass in an undetected 
co mponent (e.g., in the form of very cold molecular gas 
iPfenniger fc Reva5l2005l ). but its mass must be <C M b 
from eq. [5] so that the data stay within the scatter. Con¬ 
sequently, the stars and cold gas that we detect appear to 
compose the bulk of the baryonic mass in disk galaxies. 

4. STELLAR MASSES FROM THE BTFR 

Here we employ the BTFR as a method to estimate 
stellar mass. Rather than trust population synthesis 
models, we ask what the stellar mass must be in order 
for a galaxy to he on the BTFR. This would obviously 
be a completely circular exercise if we were to use the 
BTFR calibrated with star dominated galaxies as in the 
previous section . Fortunately , recent work on gas dom¬ 
inated galaxies (iBeeum et al. 1120081 iSwaters et al. I [20091 
iStark et aI1l200& iTrachternach et al.l l2009f i has made it 

possible to obtain an absolute calibration of the mass 
scale of the BTFR tha t is effec tively independent of stel¬ 
lar mass estimates (lMcGaughll2012lh 
The ideal situation would be a calibration of the BTFR 
with galaxies made purely of gas. In practice, this is 
impossible. Nevertheless, it is possible to select galaxies 
that are sufficiently gas rich that the error in the stellar 
mass becomes a minor contrib utor to t he er ror budget. 
T his occurs when M n > AT* (lMcGaughll2011lh 
I McGaugh] (I2Q12T) discusses fits to the BTFR for several 
distinct samples of gas rich galaxies. Combined, these are 
consistent with the BTFR 

\ogM b = 1.67+ 4 log V). (6) 

In order to estimate the stellar mass, we assume that star 
dominated galaxies exactly follow eq. [6] This seems like 


a reasonable assumption given the results of the previous 
sectio n. However, the gas rich galaxies fit by iMcGaughl 
(]2012T) are much slower rotators than the sample of ob¬ 
jects here (Fig. [6]). It is therefore a big extrapolation to 
apply eq. HU to the bright galaxies in Table Q] 

We first use the measured rotation velocity and eq. [6] 
to estimate of the total baryonic mass. The stellar mass 
then follows by subtracting the gas mass: 

M* = M b - M g . (7) 

The results of this computation are reported in Ta¬ 
ble [3] The first column of Table [3] gives the name of 
each galaxy, as in Table [D Column 2 gives the stel¬ 
lar mass as estimated by the [3.6] luminosity assuming 

Y* 3 ' 6 ^ = 0.47 Mq/Lq. This is provided for direct com¬ 
parison to the stellar mass estimated with eq. 0 which 
is given in column 3. The formal uncertainty in the stel¬ 
lar mass is estimated by propagating the uncertainty in 
Vf. distance, and the gas rich BTFR calibration. The 
corresponding mass-to-light ratios in the V, /, and [3.6] 
bands are given in columns 4-6. 

Estimating stellar masses and mass-to-light ratios with 
this procedure is very demanding on data quality. Start¬ 
ing from Vf which is typically of order ~ 100 kms -1 , 
we first raise this number to a high power. We then 
multiply by an empirically determined constant of the 
same order. This gives a large number appropriate to 
the baryonic mass of a galaxy. We then subtract another 
large number, the gas mass. This procedure invites the 
propagation of errors. The stellar masses we obtain are 
nevertheless quite reasonable (Table[3|). We should, how¬ 
ever, not be surprised if this method sometimes returns 
unphysical results, like negative stellar masses. 

Indeed, negative stellar masses are inferred four times, 
in DDO 168, D500-2, NGC 2366, and IC 2574. These 
four galaxies are the points that lie slightly above the 
line in Fig. [2] In no case is the stellar mass significantly 
negative: all four are within 1 <t of being positive. 

The opposite situation also happens. DDO 154 is in¬ 
ferred to have an unreasonably large mass-to-light ratio 
it is the lowest point in Fig. [2] However, its stellar 
mass is non-zero by less than 2 a, so its large inferred 
mass-to-light ratio is not meaningful. 

The few extremal cases are an indication of the suc¬ 
cess of the procedure. If we attempt the same exercise 
using line-widths instead of resolved rotation curves, we 
routinely find negative masses. We often also find super- 
maximal mass-to-light ratios, which do not occur in Ta¬ 
ble [3] The lower scatter in the BTFR obtained with 
resolved rotation curves is essential to the method. 

One should nevertheless treat the results for individual 
galaxies with caution. The data are statistically indica¬ 
tive, but it is very hard to have confidence in any one 
stellar mass. Given the large extrapolation from the gas 
rich regime, the stellar mass of the brightest galaxies is 
extremely sensitive to the slope of the BTFR (see H6.5I) . 

We can check the inferred mass-to-light ratios against 
the expectations of population synthesis models. This 
is done in Figure [3] where we show the mass-to-light 
ratios from Table [3] as a function of color. Figure @] shows 
histograms of the BTFR mass-to-light ratios. 

The mass-to-light ratios estimated from the BTFR are 
consistent in amplitude with that expected from popu- 
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Table 3 

Stellar Masses and Mass-to-Light Ratios from the BTFR 


Galaxy 

M*[3.6] a 

M, (BTFR) 

(10 9 M 0 ) 

y'V 

* 

I 

(M 0 /L 0 ) 

x [3.6] 

DDO 154 

0.0341 

0.248 ±0.142 

4.81 ± 2.76 

8.44 ±4.84 

3.42 ± 1.96 

D631-7 

0.0458 

0.156 ±0.147 

2.89 ±2.72 

3.48 ±3.27 

1.60 ± 1.51 

DDO 168 

0.0881 

-0.0431 ± 0.0866 

-0.26 ±0.53 


-0.23 ±0.46 

D500-2 

0.186 

-0.158 ±0.413 

-0.52 ± 1.35 

-0.70 ± 1.83 

-0.40 ±1.05 

NGC 2366 

0.336 

-0.365 ± 0.419 

-0.80 ±0.92 

-0.99 ± 1.13 

-0.51 ±0.58 

IC 2574 

1.10 

-0.373 ± 0.435 

-0.36 ±0.42 

-0.39 ± 0.46 

-0.16 ±0.19 

F563-1 

1.34 

2.36 ±2.73 

2.09 ± 2.42 

1.89 ±2.19 

0.83 ±0.95 

NGC 2976 

1.50 

2.23 ±0.53 

1.98 ± 0.47 

2.13 ±0.51 

0.70 ±0.16 

F568V-1 

1.97 

5.71 ±3.85 

2.42 ± 1.64 

2.54 ± 1.71 

1.36 ± 0.92 

NGC 1003 

2.78 

0.41 ± 1.12 

0.15 ± 0.40 


0.07 ±0.20 

F568-1 

3.33 

2.98 ±3.13 

0.96 ± 1.01 

1.00 ± 1.06 

0.42 ± 0.44 

NGC 7793 

3.55 

5.59 ± 1.43 

1.87 ±0.48 

3.10 ±0.79 

0.74 ±0.19 

UGC 128 

5.30 

4.62 ±4.58 

1.03 ± 1.02 

1.10 ± 1.09 

0.41 ± 0.41 

NGC 2403 

5.71 

11.9 ±2.8 

3.08 ± 0.72 

2.38 ±0.55 

0.98 ± 0.23 

NGC 925 

7.42 

1.42 ±2.00 

0.17 ±0.24 

0.17 ±0.24 

0.09 ± 0.12 

NGC 3198 

14.7 

21.2 ± 3.9 

1.71 ± 0.32 

1.34 ±0.25 

0.68 ± 0.13 

NGC 3621 

15.5 

6.61 ±3.86 

0.98 ±0.57 

0.93 ±0.54 

0.20 ±0.12 

NGC 3521 

42.8 

55.6 ±24.8 

3.57 ± 1.59 

2.40 ± 1.07 

0.61 ± 0.22 

NGC 3031 

49.6 

62.3 ±22.5 

2.70 ±0.97 

1.61 ±0.58 

0.59 ± 0.26 

NGC 5055 

63.5 

37.8 ± 14.8 

1.44 ±0.56 

0.94 ±0.37 

0.28 ±0.11 

NGC 2998 

75.0 

47.9 ± 13.0 

0.64 ±0.17 


0.30 ±0.08 

NGC 6674 

108. 

94.1 ±22.9 

1.49 ± 0.36 


0.41 ± 0.10 

NGC 7331 

122. 

153. ± 32. 

4.00 ± 0.84 

2.28 ±0.47 

0.59 ±0.12 

NGC 801 

126. 

42.9 ± 15.4 

0.60 ±0.21 


0.16 ± 0.06 

NGC 5533 

142. 

106. ± 24. 

1.69 ±0.38 


0.35 ± 0.08 

UGC 2885 

220. 

304. ± 69. 

1.70 ±0.38 


0.65 ± 0.15 


a Stellar masses assuming Y = 0.47 Mq/Lq are included for comparison. 


lation synthesis. This was not guaranteed: the BTFR 
calibrated by gas rich galaxies provides a constraint that 
is independent from the population models. While it was 
clear already from Fig. |T] that this would follow, it need 
not have. Indeed, it would not have if we employed a 
model much different from that of eq. 0 of which there 
are many examples in the literature. 

Nevertheless, the amplitude, slope, and scatter in T* 
expected form population synthesis are consistent with 
the gas rich BTFR. While the data are consistent with 
the CMLR determined photometrically, they are not pre¬ 
cise enough to further constrain it. That is, we cannot 
use these data to measure the color slopes of the CMLR. 

Indeed, Fig. [3] does not show a strong dependence of 
Y* on color. As expected, the dependence of Y* on B V 
becomes weaker as one goes from optical to NIR. There 
is no hint of a slope at all at [3.6]. As expected from a 
population perspective, the NIR Y* from the BTFR is 
effectively constant. 

If the dependence of Y* on B — V is weak, it is vir¬ 
tually non-existent for the other colors. There is a hint 
in the data that Y \ might become larger for very red 
V — I. Otherwise, Y* does not appear to be particu¬ 
larly sensitive to either V — I or J — K. This is not 
surprising . I t was a lrea dy c lear from the photometry 
(|McGaugh fc Schombertl 120141 1 that B — V contained 
most of the information about Y*, with V — I carry¬ 
ing only a small additional amount beyond that. There 
is no us eful infor mation abou t Y* in J — K (see also 
iScliombert fc McGaughll2014alh 

The histogram of Y* (Fig. 2]) shows scatter about a 
typical value. The widths of the histograms basically 
reflect the propagation of errors: the data are consistent 
with a narrow intrinsic scatter 1 43.51) . Nevertheless, it is 


worth noting that there is a clear and well defined peak 
at a preferred value in the NIR: Y* = 0.45 Mq/Lq. 

The top row of Fig. [Ijshows the raw distribution of Y*, 
and does not consider the expected variation with color. 
Since there is no expected variation with color in [3.6], 
the peak is clearer there than it can be in either V or /. 
In the bottom row of Fig. [4] we show the histogram of 
Y* relative to that expected from eq. [3] Here it can be 
seen that the data do peak near the value expected from 
population synthesis: there is value in the color informed 
Y* estimate of the CMLR. 

An obvious consequence is that the luminosity in any 
band is a good approximation of the stellar mass. This 
approximation can be improved with use of color infor¬ 
mation. The amount of improvement is limited by the 
intrinsic scatter in the CMLR. Both the scatter and the 
color term are minimized in the NIR. This makes the 
NIR luminosity the best available proxy for stellar mass. 

The results here are consistent with the results of 43.51 
The stellar population estimate of the BTFR gives a rela¬ 
tion indistinguishable from the separate gas rich calibra¬ 
tion. It makes little difference which method we choose. 
Either the data fall exactly on the BTFR and give the 
scatter seen in Fig. [3j or the galaxies follow exactly this 
lines of eq. [3] in Fig. [3] with the scatter seen in Fig. [I] 
There has to be some scatter in Y*. leaving precious lit¬ 
tle room for intrinsic scatter in the BTFR. 

It remains unclear why the BTFR is so well posed as 
a relation between baryonic mass and the flat rotation 
velocity. It does not follow from structure formation 
simulations without further assumptions. It has become 
common to i nvoke feedback from star formation in this 
context le.g.. iGovernato et all 120071 : lDuttonll2012th but 
it is hard to see how the chaotic processes of supernovae 
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Figure 3. Stellar mass-to-light ratios (Table [3} inferred from the gas rich calibration of the BTFR. Mass-to-light ratios in V (left column), 
I (middle column), and [3.6] (right column) are plotted as a function of color in B — V (top row), V — I (middle row), and J — K s (bottom 
row). Data with uncertainties in stellar mass < 50% are shown as solid points. Less accurate data are shown as open circles. Galaxies 
with formally negative stellar masses are shown as 2cr upper limits (grey triangles). The mass-to-light ratios predicted by eq. [3] axe shown 
as lines in the panels of the top row. 


returning energy to the ISM can result in the negligible 
scatter observed in the BTFR. The only theory to an¬ 
ticipate the specific M h-Vf relation that we observe is 
MOND riMilgromlll983h . 

Regardless of the physics underlying the BTFR, it is of 
enormous utility to have such a relation. We emphasize 
that the BTFR constructed here is empirical, making no 
reference to either dark matter or MOND. It relies only 
on observational data and stellar population synthesis 
models to inform the CMLR. 

5. A CONSTRAINT ON THE IMF 

Returning to the matter of stellar mass, there is ad¬ 
ditional information in the BTFR. The normalization of 
the BTFR calibrated by gas rich galaxies is effectively 
independent of stellar mass estimates. In contrast, the 


population synthesis mass estimates are very sensitive to 
the IMF. That the two are consistent implies that we 
are not far wrong about the IMF. Nevertheless, we can 
query the data to see what mass-to-light ratio is pre¬ 
ferred. Rather than take the population synthesis mass 
at face value and fit the BTFR, we adjust the stellar 
mass of bright galaxies to find the best agreement with 
the BTFR previously fit to gas rich galaxies. This leads 
to a constraint on the IMF fetark et alJl2009f ) . 

In order to find the ideal match with the BTFR, we 
adjust the zero point in eq.[3]while holding the color slope 
fixed. We their find t he zer o point that best matches the 
BTFR of IMcGaughl That is, we determine the 

value Aa in 


log T* = a + b(B — V) + A a 


( 8 ) 
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Figure 4. Histograms of stellar mass-to-light ratios inferred from the BTFR in V (left column), I (middle column), and [3.6] (right 
column). The top row shows the measured mass-to-light ratios while the bottom row shows T* relative to the population synthesis 
expectation of eq. QJ] (dashed lines). The bin size and plot scale are the same for all filters within each row. The data are consistent with 

t[3.6] _ q 45 Mq/Lq with modest intrinsic scatter. 


that best matches the gas rich BTFR. 

The net effect of this exercise is a small shift of the 
intercept in each band. The result is A a = —0.015 in V, 
+0.06 in /, and —0.019 in [3.6]. This is what it takes to 
make the dashed lines in the bottom row of Fig. [T] fall on 
top of the solid lines. 

These shifted zero point are not significantly differ¬ 
ent from the original population synthesis values. For 

example, the NIR value changes from tJ 3 ' 6 ^ = 0.47 to 
0.45 Mq/Lq. Similarly, in the galaxies for which we 
have Ji-band data, the mean shifts from Tf = 0.60 to 
0.57 Mq/Lq (Aa = —0.022). This corresponds to a ratio 

Tf /Ti 3 '= 1.27, consistent with the 1.29 we estimated 
photometrically (llMcGaugh fc Schombertll2C)14D . 

This nice consistency may seem mundane, but is ac¬ 
tually a rather important observation. That Aa is small 
implies that the population models are nearly correct, 
provided they adopt the right IMF. The data are consis¬ 
tent with a Kroupa or Chabrier IMF. The Salpeter IMF, 
as us ually e mploy ed, is a bit too heavy (see discussion in 
IStark et al.l[2009l) . This is consistent with what we know 
about the IMF locally, in which the numbers of low mass 
star s does not continue to increas e as a Salpeter power- 
law (lKroupalll998t [Chabriedl2003D . 

6. COMPARISON TO OTHER METHODS 

The absolute calibration of stellar mass is of obvious 
interest. In this section, we compare the normalization 
found from the BTFR to other methods in the literature. 
A largely consistent picture emerges, with one important 
exception. 


Table 3] summarizes various estimates of the stellar 
mass-to-light ratio of disk galaxy stellar populations. 
Column 1 gives a name for each method considered as 
described in column 2. Columns 3-6 give typical mass- 
to-light ratios in the V, /, K , and [3.6] bands. We in¬ 
clude the Tv-band here for comparison with [3.6]. These 
two bands are nea rly interchangeable photometrically 
(ISchombert fe McGaughlf2014bl la). We convert between 
them assuming the ratio of mass-to-light ratios is con¬ 
stant so that results limited to one or the other band can 
be directly compared. Column 8 gives the reference for 
each T* calibration. 

Our result from ]J2]is given first. These are consistent 
with our previou s attempt to employ this method 
([Stark et al.l I2009H . It is also consistent with many, 
but not all, population synthesis estimates. The NIR 
Y* are hardly different from eq. [3] The tweak in the 
H-band is also negligible, but that in the /-band is 
more noticeable. These are sensitive to color, and 
are estimated at B — V = 0.6 to provide a uniform 
comparison. We do not consider the difference between 
the BTFR T* ss 1.4 Mq/Lq and the value from eq. [3] 
of ss 1.2 Mq/Lq to be particularly significant. 
It perhaps serves more as an illustration of just how 
sensitive the method is, as this is the shift (Aa = 0.06) 
required to match the two nearly indistinguishable 
lines in the bottom middle panel of Fig. [T] Indeed, the 
positive Aa in the /-band (when the other filters are 
slightly negative) may simply be the result of missing 
data in that filter for some of the brighter galaxies. 
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Table 4 

Mass-to-Light Ratios Estimates 


Calibration 

Method 

TO 

<Ti> 

<Tf> 

(Yi 3 - 61 ) 

Ref. 

Gas Rich 

BTFR 

1.38 

1.42 

0.57 

0.45 

1 

DiskMass® 

Vertical a z 

0.74 

0.65 

0.31 

0.24 

2 

DiskMass® 

DM-adjusted 

0.57 

0.50 

0.24 

0.19 

3 

Milky Way 

Star counts 

1.5 

1.2 



4 

Milky Way 

Terminal vt 

1.53 

1.22 

0.63 

0.49 

5 

Milky Way 

Vertical K z 

1.64 

1.31 

0.66 

0.51 

6 

LMC 

Star counts 



0.65 

0.5 

7 

Boll 

Popsynth 

1.43 

1.25 

0.73 

0.62 

8 

Portinari 

Popsynth 

1.32 

1.11 

0.50 

0.41 

9 

Zibetti 

Popsynth b 

1.07 

0.76 

0.21 

0.14 

10 

Into 

Popsynth b 

1.19 

0.99 

0.41 

0.33 

11 

Eq. [3] 

Popsynth 

1.43 

1.24 

0.60 

0.47 

12 

Schombert 

Popsynth 



0.65 

0.5 

13 

Meidt 

Popsynth 



0.77 

0.6 

14 
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Note. — Color dependent population synthesis models are evalu¬ 
ated at B — V = 0.6. To compare the K s and [3.6] bands we assume 

= 1.29Y* 1 ' 6 ' llMcGaugh fe SchomberHl2014l ). with the exception of 

the BTFR method which independently suggests Y(/ = 1.27 y[, 3 ' 6 ^ (see 
text). 

® For comparison with other methods, the V and /-band values are 
approximated by scaling from the A'-band measurement with eq. [3] 
b These models do not give self-consistent stellar masses unless the 
NIR mass-to-light ratio is adjusted to larger values: Y// Rt 0.58 and 
y[, 3 6 ^ Ri 0.45 Mq/Lq (see discussion in lMcGaugh fc Schombertll2014l) . 


6.1. The DiskMass Survey 

One of the most important works constraining the 
stellar masses of disk galax ies is the DiskMass survey 
(jBershadv et all 1201(1120111 ). DiskMass observes veloc¬ 
ity dispersions across the disks of face-on galaxies. It 
then assumes a disk thickness that is statisticall y consis¬ 
tent with observations of edge-on galaxies dKregel et al.1 
120051) in order to estimate the stellar mass profile. This 
method measures the vertical restoring force to the disk, 
and is independent of the IMF. 

The results of the DiskMas s survey are summarized 
in the second and third rows. IMartinsson et all (l2013al ) 
found a near constant NIR mass-to-light ratio of (Tf-) = 
0.31 ± 0.07 Mq/Lq. The small scatter and near¬ 
constancy of the mass-to-light ratio is consistent with 
our findings here. The normalization differs by nearly a 
factor of two. 

The rather sub-maximal disks found by the DiskMass 
survey imply that the dark matter halo contributes a 
non-negligible amount to the observed velocity disper¬ 
sions. This will inevi tab ly drive the implied mass of the 
disk down further. iSwaters et al.1 (120141) make a first 
estimate of this effect, revising the mean mass-to-light 
ratio down to (T(f) = 0.24 M 0 /L 0 . This is a factor 
of ~ 2.4 times lower than what we find here. The ef¬ 
fect of the dark matter halos may be compensated some- 
what by the prese n ce of s uper-thin disks in some galaxies 
(jSchechtman-Rook fe Bershadvll2014D . We will return to 
the factor of ~ 2 discrepancy between the DiskMass and 
BTFR results after discussing other constraints. 


6.2. The Milky Way 

Our own Milky Way provides a unique environment 
in which we can hope to both count stars directly and 
constrain their mass from kinematics. This is a rich 
field that is rapidly becoming ric her; we cons ider here 
only a few relevant results. iFlvnn et al.l (|2006l ) use star 
counts from two independent surveys to constrain Galac¬ 
tic structure. Among other things, they estimate the 
mass-to-light ratio of the stellar disk of the Milky Way 
to be T)f = 1.5 ± 0.2 and T [ = 1.2 ± 0.2 M^/ Re¬ 
star ting from an initially exponential disk, iMcGaughl 
(j2008D deformed the radial mass distribution of the 
Galactic disk in order to fit the terminal veloc¬ 
ity curve in _the_ fourth quadrant ( Luna et al. 21)0 0: 
iMcClure-Griffiths fc Picket! 120071 ). The result is a mass 
model with bumps and wiggles reminiscent of those ob¬ 
served in other spirals. Indeed, the Centaurus arm 
is prominent in the inferred Milky Way mass profile. 
Matching the details of the terminal velocity curve re¬ 
quires the disk to be nearly maximal. The stellar mass 
ob tained i n this way is ~ 2% greater than estimated by 
IFlvnn et al.l (120061) : M* = 5.48 x 10 10 M 0 for R 0 = 8 
kpc. To estimate the mass-to-light rati o, we adopt 
t he to tal /-band luminosity estimated by IFlvnn et alJ 
(|2006D, Li = 4.5 x IQ 10 L 0 , and the Af-band luminosity 
from iDrimmel fe Spergell (12001D . Lk = 8.6 x 10 10 L 0 , 
which includes a somewhat uncertain correction for the 
bulge fraction. The result is Y* = 1.22 Mq/Lq and 
T f = 0.63 Mq/Lq 

iBovv fe Rixl (120131) have measured the vertical restor¬ 
ing force to the disk K z (\Z\ < 1.1 kpc) over a wide 
range of Galactocentric radii. They obtain a mass for 
the stellar disk of M* = 4.6 ± 0.3 x 10 10 M 0 . Com- 
bining this with the disk-only luminosity estimated by 
IFlvnn et al.l (I20061) . L f sk = 3.5 x IQ 10 L 0 and Lf sk = 
6.9 x 10 10 L 0 from IDrimmel fc Soergel (|2001l ) yields 
Ti = 1.31 ± 0.09 Mq/Lq and Tf = 0.66 ± 0.04 M 0 /L 0 . 
The bulge component i s exclu ded from this calculation, 
as the measurement of lBovv fc Rbd ( 2013 ) is specific to 
the Galactic disk. The uncertainty in the mass-to-light 
ratio here is o nly that from the error in stellar mass 
quoted bv IBovv fe Rixl (120131) , and does not include the 
uncertainty in the luminosity. The latter translates to 
roughly ±0.2 in T* in absolute terms. 

The three independent methods just discussed, star 
counts, the radial force in the disk, and the vertical force 
from the disk, all find mass-to-light ratios for the Milky 
Way that are consistent within the errors. They are 
also consistent with the result from the BTFR. They 
are not consistent with the low mass-to-light ratios im¬ 
plied by DiskMass, despite t he similarity of the DiskMass 
approach to that of IBovv fc Rixl (2013 ). 

6.3. The LMC 

lEskew et al.l (12012 ) perform star counts in the LMC 
with NIR data. They sum up the mass of the stars 
observed to find tJ 3 6 ^ = 0.5 M 0 /L 0 . This method 
depends some on the IMF, as one does have to extrap¬ 
olate to unseen, low mass stars. As in the Milky Way, 
the stars are counted directly, albeit to a brighter limit, 
so the assumption on the IMF is less strong than in 
population synthesis models. 
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V f (km s 1 ) V f (km s *) 

Figure 5. The BTF R constructed from the galaxies in the DiskMass survey IIBershadv et al.li2010 ). Both panels are identical along the 
abscissa, with Vf from IMartinsson et al.l ( 2013ta i. They differ in stellar mass along the ordinate. The left panel shows the baryonic mass 
for stellar masses with a constan t if-band mass-to - light r atio from population models = 0.6 Mq/Lq: McGaugh & Schombcrt 2014) 

pl us the gas masses tabulate d by Martinsson et al. (2013a). The right panel uses the same ga s masses, but adopts the stellar masses found 
by Martinsson et al. (2013a). The BTFR calibration of the gas rich galaxies (McGaugh 2012) is shown as a solid line in both panels, with 
its lcr uncertainty as dashed lines. 


6.4. Population Synthesis 

Stellar population synthesis is a well developed field. 
It is possible to use our knowledge of stellar evolution to 
construct detailed models for the composite stellar popu¬ 
lations of spiral galaxies for various assu med star forma¬ 
tion histories and metallicit ies fe.g.. iBruzual fe Charlotl 
120031: iLe Borgne et a.l.ll200dl) . For an assumed IMF, such 
models can be used to estimate the mass-to-light ratio of 
a population. 

One expects that as stars age, the population will red¬ 
den and the mass-to-light ratio will increase. Many mod¬ 
els have been constructed to quantify this behavior. A 
few relevant ones ar e noted in T able Hi 

The models of Bel l et, al.l ([200^_ ! _jPortinajiet_alJ 
(|2004fl . IZibetti et al. (j2009l ). andllnto fc Port, inaril ( 20131) 
have b een explored in detail by iMcGaug h fc Schomberti 
(120141) . These all give comparable mass-to-light ratios 
in the optical portion of the spectrum: 1 < T/ < 
1.5 Mq/Lq evaluated at B — V = 0.6. Small differences 
there appear to be attributable largely to modest differ¬ 
ences in the assumed IMF. However, the models have 
very different spectral energy distributions in t he NI R, 
w here agreement breaks down. In terms of T iBell et al.1 
(l2003h estimate T f = 0.73 M 0 /L 0 for B — V = 0.6 
while IZibetti et all (I2009T) give Tf- = 0.21 M 0 /L 0 . This 
large discrepancy persists if we corre ct for difference in 
the IMF: scali ng up the Zibetti et al.1 (120091) estimate to 


match that of lBell et al 


(120031) in the U-band has little 
impact in the NIR: Tf = 0.21 —> 0.28 M 0 /L 0 . This 
is st ill a f actor of 2.6 lower than the 0.73 M 0 /L 0 of 
iBell et all (120031) . The models do not produce c onsistent 
s pectral energy distributions dKriek et ahll201f)D . 

IMcGaugh fc Schomberti (120141 ) used multi-wavelength 
photometric data to test these various models. Requir¬ 
ing self-consistency - i.e., that the same mass of stars 
produce the observed luminosity in each band for every 
galaxy — leads to eq. [3] This stipulation obliges all mod¬ 
els to fall in the range 0.4 < tI 3 ^ < 0.5 M 0 /L 0 . 

The reasons for the difference between models are 


many, but the primary issue appears to be the treat¬ 
ment of TP-AGB stars. The luminosity of this compo- 
nent i s grossly exaggerated in the models of lZibetti et al.1 
(120091) . The TP-AGB stars, as modeled, produce a lot 
of NIR photons without contributing much in the opti¬ 
cal nor representing much stellar mass. Consequently, 
the NIR mass-to-light ratios are u nderesti ma ted. A sim¬ 
ilar conclusion was reached by IZibetti et al.l (I2013D , who 
sought but did not find the predicted spectral features of 
TP-AGB stars. 

We constructed our own population models in 
iSchombert fc McGauglil ( 2014a ). These models have an 
empirically motivated distribution of stellar metallicities 
as well as different ages corresponding to a variety of 
star formation histories. The metallicity distribution 
weakens the dependence of the NIR mass-to-light ratio 
on color relative to models built with a single metallic- 
ity. Most of model parameter space is consistent with 
T! 3 - 61 = 0.5 Mq/Lq. Independently. IMeidt et all (|2014T ) 

obtain Ti 3 6 ^ = 0.6 M 0 /L 0 albeit with a stronger color 
dependence. These values are consistent with the major¬ 
ity of methods. The exception to this is the DiskMass 
result, which we examine further in the next section. 

6.5. The DiskMass BTFR 

Examination of Table [I] suggests that nearly all meth¬ 
ods are consistent with a NIR mass-to-light ratio T^ s=s 

0.57 and T* 3 ' 6 ^ ss 0.45 M 0 /L 0 . Since the population syn¬ 
thesis result for a Kroupa IMF is consistent with the nor¬ 
malization independently required by the gas rich BTFR, 
it is tempting to conclude that the problem of stellar 
mass in disk galaxies is essentially solved. Unfortunately, 
the method of vertical velo city dispersions, as employe d 
by the DiskMass project (|Bershadv et al.l 1201(1 120111) . 
does not yield a consistent result. 

Since DiskMass provides a result that is independent 
of the IMF, it deserves more consideration than the mu¬ 
tual consistency of the various population synthesis mod¬ 
els. We can maintain that consistency while rescaling 
the IMF to match the DiskMass scale with the trans- 
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Figure 6. The BTFR from the [3.6] sample (solid red circles) compared to the gas rich gala xy sample of McGaugE] l| 2012[. open circl es') 
and the DiskMass data (squares). The recently discovered, very low mass dwarf galaxy Leo P lIRhode et al.l 20T^TGiovaneIIi et al .112013 1 is 

also shown. In the left panel, the [3.6] data are shown with the population synthesis mass-to-light ratio xi 3 ' 6 ^ = 0.47 Mq/Lq and 0.12 dex 
uncertainty. In the right panel, the DiskMass data are shown as in Fig. [5] The stellar masses of the other two data sets are scaled down 
to match the DiskMass scale (y[?’ 6 ^ = 0.24 Mq/Lq for the [3.6] data). This makes little difference to the gas rich samp le , but the [3.6] 
sample shifts to match the DiskMass data. The solid line in both panels is the fit to the gas rich data alone from iMcGaiigiit <2012 1. A fit to 
the combined gas rich plus [3.6] sample is shown as the dashed line in the left panel. In the right panel, the dashed line is a fit to the gas 
rich plus DiskMass sample. A tiny adjustment in slope (see Table l5|l over the six decades in mass illustrated here suffices to accommodate 
the DiskMass data with the gas rich data. 


lation A a = —0.29. So conceivably we simply need to 
adopt a different IMF in the population modeling. While 
this may be acceptable form a population synthesis per¬ 
spective, it is harder to reconcile with Milky Way con¬ 
straints. We check here whether we can reconcile it with 
the BTFR. 

We construct the BT FR from the DiskMass data 
(|Martinsson et al.l 12013 b. al for two choices of mass- 
to-light ratio. In one case, we use t hat mea¬ 
sured by DiskMass ijMartinsson et all l2013af) . In the 
other, we adopt a constant mass-to-light ratio = 
0.6 Mq/Lq as suggeste d by po pulation synthesis models 
(jMcGaugh fe Schombertll2014D . The result is shown in 
Fig. [5] Unsurprisingly, the two BTFR are shifted with 
respect to each other due to the difference in the adopted 
mass-to-light ratios. 

The BTFR constructed using the population synthesis 
value Tj- = 0.6 Mq/Lq places the DiskMass galaxies di- 
re ctly on the BTFR fit independently to gas rich galaxies 
bv iMcGauehl (120121) . This is not guaranteed to happen. 
The gas rich galaxies are, by and large, much slower rota¬ 
tors than the bright spirals of the DiskMass sample. The 
extrapolation of that fit will not work if the mass-to-light 
ratio is discordant. 

Indeed, the data for the mass-to-light ratios deter¬ 
mined by DiskMass are systematically shifted off the gas 
rich BTFR. Most of the data parallel it, falling ~ ler be¬ 


low the fit line. The DiskMass data are consistent with 
a slope 4 BTFR, but with a lower normalization. One 
possibility is thus that the DiskMass Y* are too low due 
to some unknown systematic error, and simply need to 
be shifted upwards to be consistent with other results. 

Looking closely at the DiskMass BTFR (right panel 
of Fig. O, one may note that the data almost bifurcate 
into two parallel relations. Most of the data fall low, 
but a half dozen galaxies stand to one side, consistent 
with the gas rich BTFR. We do not know why the data 
might split like this, and suspect it may be indicative of 
a systematic error. This bifurcation disappears when a 
higher, constant IL-band mass-to-light ratio is adopted. 
The scatter also decreases, but this is an inevitable result 
of setting the inclinations of the DiskM ass galaxies wi th 
the iL -band LTFR in the first place (fMartinsso n et all 
I20l3al) . As Y* grows, the gas becomes less important to 
Mb, and one eventually just recovers a scaled version of 
the LTFR. 

So far, we have kept distinct the various samples: gas 
rich galaxies, the [3.6] sample of this paper, and the 
DiskMass data. Another approach is to make a fit of the 
BTFR with the combined data. The [3.6] and DiskMass 
data are consistent with each other for the same choice of 
mass-to-light ratio, so we make two combined fits. In one, 
we simply combine the [3.6] data as presented here with 
the gas rich data from iMcGaughl (120121) . In the other, 
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Table 5 

BTFR of Combined Samples 


Sample 

X 

log A 

3.6 + gas rich 

4.04 ±0.09 

1.61 ±0.18 

DiskMass + gas 

3.87 ±0.10 

1.80 ±0.21 


Note. — Fits of the form log Mb = log A+ 
x log Vf. 


we combine the DiskMass data with the gas rich data, 
scaling the stellar masses of the gas rich galaxies down 
by a factor of two to be consistent with the DiskMass 
scale. These fits are shown in Fig. [6] with fit parameters 
given in Table [5] 

In combining the [3.6] and gas rich samples, we have 
been careful t o exclude du plica te galaxies. The gas rich 
galaxies from IMcGaugh (l201 2i) are largely independent 
of the sample here, but th ere is some ov erlap. Of the 
47 galaxies considered in iMcGauehl (l2012h and the 26 in 
TablelH there are 7 objects in common: D631-7, DDO 
168, IC 2574, F568-V1, F568-1, NGC 2403, and NGC 
3198. These seven are plotted (and fit) only once in 
Fig. El as part of the [3.6] sample. 

The combined fit obtained in this fashion is consistent 
with the previous fit to the gas rich data alone. It is 
slightly steeper, as we have retained the stellar popula¬ 
tion estimate of xi 3 ’^ = 0.47 Mq/Lq, which is ever so 
slightly higher than the 0.45 of a straight extrapolation 
from the gas rich BTFR. The fit lines are hardly distin¬ 
guishable. 

Indeed, the [3.6] data follow the extrapolation of the fit 
to the gas rich galaxies remarkably well. This happens 
with zero effort: the location of the galaxies from Table [l] 
in the BTFR plane is determined entirely by the data and 
the mass-to-light ratio provided by population modeling. 
It is completely independent of the gas rich calibration. 

The gas rich BTFR calibration also extrapolates well 
downwards to lower masses. This holds for both star 
dominated dwarf Spheroidals (IMcGaugh fc Woll I2010H 
and for the recently d iscovered gas rich dwarf Leo P 
(iGiovanelli et al.ll2013il. In Figure [6] we adopt the dis¬ 
tance of Mc Quinn et al.l (120131! when plotting the stellar 
(|Rhode et al.ll2013H and gas (IGiovanelli et al.ll2013f) mass 
of Leo P so that the datum is independent of the fitted 
relation. Leo P falls along the extrapolation of the BTFR 
to lower rotation velocity just as the [3.6] data fall along 
the extrapolation to higher Vf. The BTFR appears to 
be well described as a single power law over nearly six 
decades in mass. 

In the combined DiskMass plus gas rich BTFR fit, the 
relation again appears to be well described as a single 
power law. The chief difference is a slightly lower slope. 
Instead of x « 4.0, we now find x ~ 3.9. 

This change of slope is driven by the lower mass- 
to-light ratios favored by DiskMass. The lower stellar 
masses adopted for the gas rich galaxies themselves make 
no perceptible difference. The baryonic masses of these 
objects were already dominated by gas; for the new as¬ 
sumption this is slightly more true. However, the gas 
rich galaxies are predominantly slow rotators, while the 
DiskMass sample is composed entirely of fast rotators. In 
order to fit the two data sets simultaneously, the BTFR 


simply pivots around the low mass end, with the fitting 
routine selecting the slope necessary to hit the data at 
the high mass end. The two lines illustrated in the right 
panel of Fig. [6] are indistinguishable over the range con¬ 
strained by the gas rich data, but the slight difference in 
their slopes amounts to a factor of ~ 2 in T* for the fast 
rotators. 

It therefore appears that a factor of ~ 2 systematic un¬ 
certainty in stellar mass persists. The natural conclusion 
of our work is that the population synthesis models are 
essentially correct for a reasonable (e.g., Kroupa) IMF. 
However, we cannot exclude the lower mass scale indi¬ 
cated by the DiskMass data. The BTFR fit to DiskMass 
is within the uncertainties of the original gas rich BTFR 
fit. 

The situation could be improved with better data. Di¬ 
rect distance measurements of many gas rich galaxies 
would considerably improve the calibration. If the small 
intrinsic scatter of the BTFR persists, then the scatter 
will come down as the data improve. The gas rich BTFR 
provides a promising new method to constrain the abso¬ 
lute stellar mass scale. 

7. CONCLUSIONS 

We have weighed stellar disks in two distinct ways. We 
have used a traditional population synthesis approach 
with u pdated color-mass-to -light ratio relations (CMLR 
IMcGaugh fc Schombertl [20141 1. Separately, we have in¬ 
ferred the stellar mass required for bright galaxies to 
lie along the baryonic Tully-Fisher relation (BTFR) cal¬ 
ibrated by gas rich galaxies. These independent ap¬ 
proaches yield consistent results. 

We construct Tully-Fisher relations in the V, /, and 
[3.6] bands together with measurements of the outer, flat 
rotation velocity from well resolved 21 cm data cubes. 
Each band gives a distinct luminous Tully-Fisher relation 
(LTFR). Each LTFR has a different slope, intercept, and 
scatter. 

The CMLR o btain ed by requiring photometric self- 
consistency by IMcGaugh fc Schombertl ( 201 41 success¬ 
fully converts the different LTFR into a consistent stel¬ 
lar mass Tully-Fisher relation (STFR). This provides 
kinematic confirmation of the ph otometric re sults. The 
slopes of the CMLR obtained by IMcGaugh fe Schombertl 
(]2014f) must be approximately correct, or self-consistent 
STFR would not follow. 

We use the scatter about the STFR to estimate the in¬ 
trinsic scatter in the CMLR. The scatter in stellar mass is 
0.16 dex in V , 0.13 in /, and 0.12 in [3.6]. The trend for 
smaller scatter towards redder wavelengths follows the 
expectation of stellar population models. Indeed, the 
modest amount of scatter in the nea r -infra red (NI R) is 
consistent with the m odels of lBell fc de .Tone] (|200lH and 
iPortinari et al.l (I2004H . who estimate 0.1 and 0.15 dex in 
the A'-band, respectively. This consistency implies that 
the variation in star formation history from galaxy to 
galaxy suffices to explain the observed scatter, leaving 
little room for other sources of scatter. This in turn im¬ 
plies that the IMF is effectively universal, at least when 
averaged over entire galaxies. 

The stellar masses estimated with our CMLR also lead 
to nice convergence in the BTFR. The BTFR obtained 
in this fashion is consistent with that o btain ed indepen¬ 
dently from gas rich galaxies (|McGaugljl2012f) . The slope 
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is steep, being very close to x = 4 in Mb oc V'j. There is 
no indication of a bend that transitions between a steeper 
slope at low mass and a shallower slo pe at high mass 
as su ggested by some models (e.g.. iTruiillo-Gomez et al.1 

l20H . 

The intrinsic scatter of the BTFR is very small, consis¬ 
tent with zero. Effectively all of the observed scatter can 
be attributed to observational uncertainties and the scat¬ 
ter in stellar mass-to-light ratios. Ignoring these, the up¬ 
per limit on intrinsic scatter is that observed at [3.6]: 0.13 
dex. This is difficult to understand in the context of dark 
matter models, in which there should be a considerab le 
amount of scatter from halo to halo (lEise nstein fc Loebl 
119961: iMcGaugh fc de Blokll 19981 [Bullock et al.ll200lH ~ 

The BTFR calibrated by gas dominated galaxies pro¬ 
vides an independent constraint on stellar mass. This 
approach also yields a typical 3.6/nn mass-to-light ra¬ 
tio of Tl 3 = 0.45 Mq/Lq, with a corresponding value 
in the A' s -band of Tf- = 0.57 M Q /L Q . The scat¬ 
ter about these typical values is 0.12 dex, consistent 
with the anticipated variation in star formation histo¬ 
ries dBell fc de .Tonal 120011 : iPortinari et all 12004 iKelsonl 

l20Ti. 

An important consequence of this result is that the 
NIR luminosity of a galaxy is a good proxy for its stellar 
mass. Workers wishing to estimate the the stellar mass 
of galaxies would do well to measure this quantity and 
adopt a constant NIR mass-to-light ratio as calibrated 
here. If a NIR luminosity is not available, an optical 
luminosity is also a good proxy, albeit with more scat¬ 
ter and some color dependence. Contrary to our initial 
hope, there appears to be little value added in fitting the 
entire spectral energy distribution of galaxies, at least 
so far as their mass-to-light ratio is concerned: most of 
the information about the mass-to-light ratio is carried 
by a single color like B — V. Further color-based correc¬ 
tions are more likely to reveal systematic uncertainties in 
stellar population models than they are to improve the 
estimate of stellar mass, at least at this juncture. 

The BTFR provides a new method to calibrate the 
absolute scale of stellar mass in rotating galaxies. It 
provides a useful constraint that is consistent with stel¬ 
lar population synthesis models. The accuracy of this 
method could be considerably improved with straightfor¬ 
ward observations, such as direct distance measurements 
of gas rich galaxies. 
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